Effects of counterion fluctuations in a polyelectrolyte brush 
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We investigate the effect of counterion fluctuations in a single polyelectrolyte brush in the absence 
of added salt by systematically expanding the counterion free energy about Poisson-Boltzmann mean 
field theory. We find that for strongly charged brushes, there is a collapse regime in which the brush 
height decreases with increasing charge on the polyelectrolyte chains. The transition to this collapsed 
regime is similar to the liquid-gas transition, which has a first-order line terminating at a critical 
point. We find that for monovalent counterions the transition is discontinuous in theta solvent, 
while for multivalent counterions the transition is generally continuous. For collapsed brushes, the 
brush height is not independent of grafting density as it is for osmotic brushes, but scales linear 
with it. 

PACS numbers: 61.25.Hq, 61.20.Qg 



I. INTRODUCTION 

Recently, correlation effects in electrostatics of highly 
charged macroions in aqueous solution have received a 
great deal of attention (see Refs. for recent re- 

views). These effects give rise to interesting phenom- 
ena that are not contained in the mean-field Poisson- 
Boltzmann (PB) theory 0]^such as charge inversion 
and like-charge attraction 2,"^. One plausible origin of 
the attraction comes from fluctuations of the counterion 
density that become correlated leading to a long-range 
attraction that is similar to van der Waals forces 0,0)l3- 
Counterion fluctuations might have important effects on 
the conformation of highly-charged polyelectrolytes 
In this paper, we study the effects of counterion fluctu- 
ations in a polyelectrolyte brush, where these effects are 
thus expected to be important. 

A polyelectrolyte brush consists of a high density of 
charged polymers with one of their ends grafted to a 
surface (see Fig^. Because of their technological ap- 
plications, such as enhanced stabilization against col- 
loidal flocculation, they have been extensively studied in 
the past both theoretically IS [H [H ElEl El 
and experimentally [itI Il9l l2Clj |. Pincus |12j| stud- 
ied polyelectrolyte brushes, based on the Alexander- 
deGennes approximation \2^ . where monomer density 
is assumed to be constant up to the brush height, and 
mean-field PB theory for the electrostatic interaction be- 
tween the charged monomers and the counterions. The 
brush height, a quantity that is accessible to experiments, 
is then determined self-consistently by minimizing the to- 
tal free energy of the brush. In a theta solvent, it sim- 
ply consists of the elastic energy arising from stretching 
the chains and the electrostatic (mean-field) free energy. 
In this picture, there are basically two regimes - weak- 
charging and strong-charging (osmotic) regimes. In the 
weak-charging regimes, the counterions are diffuse and 
the brush height scales with the degree of polymerization 
of the charged polymer TV, like h ~ N^. In the osmotic 
regime, most of the counterions are trapped inside the 



brush, whose osmotic pressure (which tends to increase 
the brush height) is balanced by elasticity of the chains 
(which tends to decrease the brush height). In this limit, 
the chains are stretched and the brush height, h ^ N and 
is independent of the grafting density. In later studies, 
this mean-field electrostatic picture has been extended 
to poor solvent, where a first-order collapse transition 
was predicted jlSj , to the quasi-neutral regime where ex- 
cluded volume effects are important |Q| , and to the self- 
consistent field theory which goes beyond the Alexander- 
deGennes approximation [l5j |. 

Early experiments on moderately charged polyelec- 
trolyte brush show that the above theory is reasonably 
valid 0, 0, 0, |23|. However, recent experiments on 
highly charged brushes have shown behavior that cannot 
be described within mean-field PB theory, especially in 
the presence of multivalent counterions. Experiments on 
the normal forces between two brushes ^] have seen de- 
creases of brush height with multivalent salt and conden- 
sation beyond that predicted by Manning theory. There 
is also evidence of attractive interactions between the 
brushes [2^. In addition, Bendejacq et al. have ob- 
served a regime in highly charged cylindrical brushes of 
annealed polyelectrolytes characterized by a continuous 
decrease in brush height with increasing brush charge. 
The nonmonotonic behavior of the brush height occurs 
even with monovalent counterions and there exists a 
novel scaling relation between the brush height and the 
concentration of added salt. Finally, a recent simula- 
tion [2^ also finds enhanced condensation and a collapsed 
brush regime even with monovalent counterions. It also 
shows for large values of grafting density, the brush height 
does not follow any of the scaling laws above. 

These observations strongly suggest that counterion 
correlation effects are important in highly charged 
brushes, especially in the presence of multivalent ions. 
In Ref. '26i|, Csajka et al. developed a scaling theory de- 
scribing strongly charged brushes using the so-called box 
model. In this model, the chains in the brush are assumed 
to be equally stretched and the counterions are assumed 
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to be distributed with constant density. In addition to 
the mean-field electrostatic free energy, a counterion fluc- 
tuation energy, estimated using the Debye-Hiickel free 
energy of electrolytes '27^, is added to reflect correla- 
tion effects. This new term tends to decrease the brush 
height. This model predicts a collapsed regime in which 
the Debye-Hiickel fluctuation pressure is balanced by the 
second virial pressure. The transition between the os- 
motic brush and the collapsed brush is predicted to be 
first order. 

In this paper, we develop a theoretical framework that 
goes beyond scaling theory and systematically study the 
role of counterion density fluctuations in the absence of 
added salt. This is done by employing a field theory for 
the counterion contribution to the free energy and sys- 
tematically expanding about the mean field counterion 
density to second order in the fluctuations 0, 113 ■ The 
chains are treated within the Alexander-deGennes ap- 
proximation of equal stretching and constant monomer 
density. We find that for strongly charged brushes, 
there is a collapse transition in which the brush height 
decreases with increasing charge on the polyelectrolyte 
chains. Unlike the scaling theory we find that the 
transition is similar to the liquid-gas transition with a line 
of first-order phase transition terminating at a critical 
point. Surprisingly, the valence of the counterions plays 
an important role: the transition to the collapse regime 
is discontinuous for monovalent counterions, while it is 
continuous for multivalent counterions. Furthermore, for 
collapsed brushes, we find that the brush height is no 
longer independent of the grafting density, as in osmotic 
brushes, but scales linearly with it. Our results are ex- 
pected to be relevant to ongoing experimental studies in 
polyelectrolyte brushes and our framework can be easily 
generalized to study the important case of the interaction 
between two highly-charged brushes.. 

The organization of the paper is as follows: In sec- 
tionHH we review the expansion about the PB mean field 
theory in order to establish notation. In section UTTl we 
specialize to the case of a single polyelectrolyte brush 
and calculate the phase diagram in both theta solvent 








and good solvent. Finally, in section HVI we discuss the 
meaning and limitations of our results. The detailed cal- 
culation of the fluctuation effects have been left to the 
Appendix. 

II. COUNTERION FREE ENERGY 

In this section, we briefly review the theoretical 
method that allows us to compute in a rigorous fash- 
ion the total electrostatic free energy that includes fluc- 
tuation contributions from the counterions. The details 
can be found elsewhere |2^ l^^. The idea is to map 
the partition function of a system of counterions inter- 
acting among themselves and with a fixed distribution 
of charges to a field theory characterized by an effective 
Hamiltonian. We can then systematically expand about 
its saddle-point solution, which corresponds to nothing 
but the mean-field PB approximation, and calculate the 
fiuctuation contributions to the lowest order. While the 
discussion here is general, we will apply this formulation 
to our model of polyelectrolyte brush in Sec. IIIII 

Consider the primitive model where a system of Nc 
point-like charged particles of charge —Ze (where Z is the 
valence) and external fixed charges with charge density 
tT(x) = en/(x), immersed in an aqueous solution with 
dielectric constant e. The electrostatic energy for this 
charged system is 

Na No 

(3E^^ = ZHb - </.(x,), (1) 

3>k ■> ] = 1 

where i^(x) — ZIb J d"^x' is the "external" electric 

potential from the fixed charge distribution and Ib = 
^i^^rp ~ TA is the Bjerrum length in water at room 
temperature, fc^ is the Boltzmann constant, and T is the 
temperature. The partition function is 

ZnA4>] = 1^ n / ^ ^-PEn.). (2) 

^' i=l •' 

where a is the molecular size of the counterions and 
f3 = l/{kBT). The partition function of Eq. Q can 
be mapped into a field theory, by switching to the grand 
canonical ensemble. We introduce a chemical potential fi 
(in units of kBT) and express the grand partition func- 
tion as 

C30 

e^'^^^^eM- (3) 

After a Hubbard-Stratonovich transformation, we obtain 
the functional representation for the system of charges in 
an external field: 



FIG. 1: A schematic picture of a polyelectrolyte brush. 



(4) 
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where is — ■^hIbZ'^, — Po^b, Po — e^'/a^, and A/q"^ = 
det[— V^] is the normahzation factor. The saddle point, 



(5) 



ss 



V'o(x), given by 
equation 



0, is determined by the 



V^[iM^)] = K2e*''"'(^)+"*(^\ (6) 

which after introducing the mean-field potential, (p{'x.) = 
— ■j'(/'o(x) — 0(x), becomes 



VV(x) 



(7) 



which is the Poisson-Boltzmann equation. This equation 
will be examined more closely in Sec. lIII^ for our model 
of polyelectrolyte brush to be introduced in Sec. Illll 

To obtain the free energy for the counterions at the 
mean-field level, we note that it is related to the Gibbs 
potential ro[(/)] = S[^q, cj)] by a Legendre transformation: 



/3Fo(n) = ro[0]-f ^ j d3xpo(x). 



(8) 



where /9o(x) is the mean- field free counterion density 
given by 



Po(x) = poe 



ii/)o(x)+0(x) 



(9) 



To capture correlation effects, we expand the action 
0] about the saddle-point, "00 (x), to second order 
in A0(x) = V'(x) - ?/'o(x): 



V] = S[(f),ipo] 



(10) 



- / / dVAVXx)K(x,y)A^(y) 



where the differential operator 



K(x,y) = 



_y2 _|_ ^2giVo(x)-|-0(x) 



(5(x 



(11) 



is the second variation of the action (j)]. Notice that 
the linear term in A?/'(x) does not contribute to the ex- 
pansion since V'o(x) satisfies the saddle-point equation 
Eq. Performing the Gaussian integrals in Eq. Q), 
we obtain an expression for the change in the free energy 
due to fluctuations of the counterions: 



total electrostatic free energy of the system can be ex- 
pressed as 



(3F,i = /?Fo + i In dct K - i In det[- V^] . 



(13) 



Thus, the total electrostatic free energy is comprised of 
a mean-field free energy and a fluctuating free energy. 
Note that the expansion parameter of the total free en- 
ergy is proportional n Hb , where l/n roughly corresponds 
to the average distance between charges. Thus, for highly 
charged systems, the fluctuating term, linear in k£b, may 
become important, and the mean- field PB free energy no 
longer provides a reasonable approximation. The com- 
putation of the fluctuating free energy relies on diag- 
onalizing the operator K(x, y), which depends on the 
mean-field solution. For counterions in a uniform back- 
ground, this term gives the usual Debye-Hiickcl free en- 
ergy, AF/V = -kBTK^/iUn) H^]. However, it might 
be difficult to calculate this term in general, and exact 
evaluation is known only for planar systems 7, 28]. In- 
deed, as we will see below, an approximation must be 
introduced in order to compute this term analytically for 
our model of polyelectrolyte brush. 



III. SINGLE POLYELECTROLYTE BRUSH 

We consider a polyelectrolyte brush (see Fig. ^1, 
formed by a long intrinsically flexible charged polymers 
grafted at one end onto an impermeable planar surface 
and immersed in a salt-free solution. We consider the 
simple case of monodisperse polymers of degree of poly- 
merization, N, of which a fraction / is charged, with 
grafting density p. We avoid the complicated issue of 
how electrostatics induces stiffness of the chains, and sim- 
ply assume that the polymers are Gaussian chains, with 
Kuhn length a which is the same molecular size of the 
counterions. 

In addition to the electrostatic free energy, we have to 
model the energetics of the polyelectrolyte brush, aris- 
ing from elasticity and its interaction with the solvent. 
Within the Alexander-deGennes approximation j23|, in 
which monomer density are assumed to be constant in 
a region of height h, all the chains are stretched by an 
equal amount and store an elastic energy per unit area 
of 



PFg/A = 



2Na? 



P. 



(14) 



where A is the area of the surface. The solvent quality is 
described in terms of the Flory-Huggins free energy per 
unit area l30l 



/3AJ'= ilndetK- ilndet[-V^], (12) 



(3Fv/A = -a^c^h + -a^c^h. 
2 6 



(15) 



where the second term comes from the normalization fac- 
tor Aq. Thus, combining Eqs. (|HJl and (|12|l together, the 



where c = Np/h is the average monomer density, v is 
the dimensionless excluded-volume parameter, which is 
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positive for good solvents and negative for poor solvent, 
and w is the third virial coefficient which is typically 
positive and of order unity. In the next subsection, we 
discuss the electrostatic free energy within our model of 
polyelectrolyte model, as outlined in Sec.|nJ In Sec. lIIIBl 
we present the phase diagram of our model. 



A. Mean-field theory and fluctuation contribution 

The fixed charge distribution for an Alexander-de 
Gennes polyelectrolyte brush is ?^/(x) — rip 8(2 — ft.), 
where uq = fc = Nfp/h and / is the fraction of 
monomers that are charged. The function Q{x) = 1 — 
6(x), where &{x) is the Heaviside step function, which 
equals 1 when a; > and when a; < 0. Due to transla- 
tion invariance in the xy-plaxie, the PB equation, Eq. ^ 
can be written as 



-—Q{z 



h). 



(16) 



Outside the brush {z > h), the right-hand side of Eq. Hlt)|) 
is zero and its solution is the same as that of a charged 
plane held at a constant potential: 



ip>{z) = (fh + 21n 



1 



-h 



A 



z > h. 



(17) 



where Lph, which will be determined below, is the poten- 
tial at the brush height and A is the Gouy-Chapmann 
length, determined by A = e'^''/^/K, which is a mea- 
sure of the extent of how diffuse the counterion distribu- 
tion is outside of the brush. 
Inside the brush we must solve 



dz^ - Z ■ 



< h. 



(18) 



Multiplying this equation with ip^, integrating once, and 
imposing the boundary conditions that the potential and 
the electric field at the grafting surface be zero, "the con- 



stant of motion" is found to be 



1 / d(p< 
2\ dz 



[e 



1] 



and we have 



(19) 



Using the continuity of the electric field and potential 
across the boundary at z — h, we find that is related 
to (fih by 



2^ 



(20) 



where we have defined ^ = 2Z / {iBn^h). Note that this 
length is independent of h. The electrostatic potential 
inside the brush, 0<(2:), is determined implicitly by 



dip' 



le-f' - 1 



'■Ph 



= = A[^<(z),(^ft] = 2^/ 1^2:. 



(21) 



In particular, the potential at the brush height, iph is 
determined by the relation 



(22) 



In the weak charging limit, it gives iph ~ h/^ <^ 1, and 
in the strong charging limit, iph ~ 1. Fig. [3 plots the 
counterion density. It is clear in the strong charging limit 
that while most of the counterions is in the brush, there 
are still some outside. This point will be crucial when we 
discuss the fluctuation contributions below. 

The mean-field electrostatic free energy per unit area 
can be calculated straightforwardly using the solution de- 
scribed above and Eq. ©. The result is 



Pfo = 



Nfp ( 2ipha 



Z 



In 



£B^h 



2^ 



(23) 



Vh 



dip 



1 + ^ 



where we have made used of the fact that the chem- 
ical potential is given hy fi = \n[ipha^ /{2£,'^£b)]- In 
the mean-field picture, electrostatics contributes an out- 
ward pressure, EmfI^) = —dfo/dh, which scales as 
IImfCi) ~ ^/i^-B^h), in the strong charging limit, and 
^Mpih) ^ 1/(^bC^) in the weak charging limit. Note that 
this mean-field pressure is always repulsive. Balance this 
electrostatic pressure against that arising from the elas- 
ticity of Eq. H14|) . we have h « £,/ Pel, for (3^1 <C 1 (strong 
charging) and h « £,/0li, for Pel ^ 1 (weak charging), 
where we have defined an important dimensionless pa- 
rameter Pel = Z'^^'^il {N f'^/'^a)^ which characterizes the 
importance of the elastic term compared to electrostatic 
pressure. These scahng regimes (as plotted in Fig. I^J 
are well-known as discussed in Ref. Here, we study 



outside brush 




20 40 



80 100 120 140 



FIG. 2: A plot of the mean-field counterion density in the 
osmotic brush limit. It is clear that most of the counterion is 
trapped inside the brush. The dashed line represents the spa- 
tially averaged counterion density inside the brush [See, Eq. 
I27I1 I. This approximation is used to compute the fluctuation 
free energy. 
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deviations from these scaling laws when counterion fluc- 
tuations are taken into account. 

To evaluate the fluctuation contribution to the free en- 
ergy, we start with the expression for the change in the 
free energy due to fluctuations of the counterions, Eq. 

m 

I3AJ^ = i In det K - i In det [- V^] , 

where the operator K for our model of polyelectrolyte 
brush is 



K(x,y) = 



2Q{z-h) 



{z-h + Xf 



5(x-y).(24) 



The derivative of AjF with respect to h can be straight- 
forwardly calculated by making use of the exact identity 
(5 In det X = TrX-i(5X: 



dh 



= 2^/ ^'^^('^'^) (25) 

2e(z - h) 



dh 



K^e~'^<Q{z-h) + 



{z-h + \Y_ 



where G(x, x) is the Green's function (inverse) of K(x, y) 
satisfying 



-^B(5(x-x'). (26) 



Unfortunately, G'(x, x') cannot be evaluated analytically 
because the counterion density inside the brush is not an 
explicit function of z [see Eq. (I21|l ]. In order to make 
analytical progress, we employ the reasonable approxi- 
mation that the counterion density inside the brush is a 
constant averaged over the brush, given by 



h Jo \ XJ 



(27) 




FIG. 3: The height versus beta for values of F = 0, 10* 
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FIG. 4: The fluctuation contribution to the free energy, Ilo{x) 
(solid line), where x = /i/^ is the rescaled height variable. The 
scaling prediction for the fluctuation pressure from equation 
I29II (dashed line) diverges as small x. 



This approximation entails neglecting the spatial varia- 
tions inside the brush. Judging from Fig. 12 this approx- 
imation should not be in serious error. 

Within this approximation, G(x, x') can be evaluated 
analytically (see Appendix), and we find 



1_ dl3/\T 
A dh 



A2 
h 



Iiih) 

1 dX 
'yh'dh 



I^{h)-I^{h) 



h^ 



2l2(/i)],(28) 



where the dimensionless integrals Xi(/i) are defined in the 
Appendix. This is the main result of this paper. It says 
that counterion fluctuations contribute an inward (nega- 
tive) pressure n/(/i) = - ^"^ 9AJ^/9/i = F^-^ no(V^), 
as plotted in Fig. ^ First, the fact that it is negative 
is consistent with the notion that the fluctuations must 
lower the free energy. The strength of this pressure is 
controlled by F = Ib/S,, which roughly scales as [Isf pY ■ 
Thus, this term becomes important for highly charged 
polyelectrolytes, as expected. 

Although the expression Eq. H28I) looks rather com- 
plicated, we can understand its meaning by expanding 
no(/i/^) in the strong charging limit: 

no(VO = -«i {^/hf^ i^/hf 

- as ii/hf^ + ■■■, (29) 

where ai ^ 0.0375 and a2 ~ 0.00833 are constants. The 
first term is the fluctuation contribution to the free en- 
ergy of a coulomb gas of counterions distributed uni- 
formly throughout the brush volume. To see this, re- 
call that in the strong charging limit most of the coun- 
terions is within the brush, with a 3D concentration 
of c '--^ no/h. This implies that the inverse of the 3D 
"screening" length Ks ^ ^JisnQ/h. Using the Debye- 
Hiickel free energy (per unit volume) /3A/ ~ — k^, we 
have n/(/i) ~ 9(A/ • h) /dh ~ h'^/^, which is the first 
term in Eq. H29|) . In fact, similar argument has been 
used by Csajka et al. 26] to derive the scaling laws. The 
second term in Eq. (|29|) may be interpret as an effective 
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second virial term arising from the coupling between fluc- 
tuations of the counterions inside the brush and outside. 
Note that this term is not included in the scaling theory 
for collapsed brushes. It is interesting to note that the 
divergence at small h of the asymptotic fluctuation pres- 
sure is purely an artifact of the approximation, and it is 
not present in Eq. H28|l . It is likely then that scaling laws 
derived for the collapsed state of a polyelectrolyte brush 
using similar approximate expressions may not capture 
the full behavior of the brush. 



B. Phase Diagram 

Now, we are in the position to discuss the equilibrium 
phase diagram of the polyelectrolyte brush. The impor- 
tant quantity, the equilibrium height of the brush, is de- 
termined by minimizing the total free energy of the poly- 
electrolyte brush, F = Fq + Fy + F^i, with respect to 
the brush height. It is more convenient to employ the 
following dimensionless variables: 



Npa^ 



(30) 
(31) 



in addition to a; = h/^, T = £b/^, and Pel = 
Z^/^£^/(N f^/^a) already deflned in the previous section. 
The variable v is a measure of the contribution to the 
pressure from the second virial coefRcient compared to 
the mean fleld pressure. Note that the mean-fleld pres- 
sure and the fluctuation pressure depends on x only. 

In Fig. [51 we display the typical phase diagram of a 
polyelectrolyte brush in contact with a theta solvent as 
a function of Pel and T. We flnd a first order collapse 
transition line, from the osmotic brush to the collapse 
brush, represented by the solid line in Fig. This line 
terminates at a critical point, which is well within the 
strong charging regime. Pel < 1- Thus, it is similar to the 
liquid-gas transition. This behavior resembles to that of 
a polyelectrolyte brush in a poor solvent, where attrac- 
tive interaction between monomers mediated by the sol- 
vent acts against the electrostatic repulsion to collapse 
the brush Here, fluctuations induce an attractive 

interaction which collapse the brush. This effect on the 
brush height can be clearly seen in Fig.|3| It is interesting 
to note that the valence of the counterions plays a crucial 
role. For typical brushes, increasing the charge fraction 
while holding other parameters fixed (as shown with dot- 
ted lines in Fig.lSJ, only the monovalent counterion goes 
through the first-order collapse. For divalent and triva- 
Icnt counterions, on the other hand, we observe a smooth 
continuous decrease in h, as shown in Fig. 1^1 where we 
plot the brush height as a function of the charge fraction 
/. As expected, correlation effects are more important 
for higher valence Z of the counterions. Even for a mod- 
erately charge fraction /, the brush height already starts 
to decrease for for divalent and trivalent counterions. It 
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FIG. 5: Typical phase diagram of a single brush in theta 
solvent (in this case. A*' — 200, a; = 1). The solid line repre- 
sents a first order phase transition between the osmotic brush 
and the collapsed brush regimes. The Pincus brush regime 
occurs for ~ 1 and is not shown. The dashed curves are 
generated holding all parameters constant except the brush 
charge fraction (arrows indicate increasing charge fraction) 
for monovalent, divalent, and trivalent counterions. 
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FIG. 6: Brush height for monovalent (solid), divalent 
(dashed), and trivalent (dotted) counterions in a theta sol- 
vent (see figure The monovalent salt height obeys the 
power law h oc J^'^", though it shows minor deviations from 
power at very strong charging (with an exponent 0.45). 



is also interesting to observe that for monovalent coun- 
terion (before the collapse transition), the brush height 
obeys the osmotic brush scaling law oi h ^ /^/^ only for 
moderate /, which crosses over to h ^ j0A5^ Deep in- 
side the collapsed regime, the brush height is determined 
by the fiuctuation pressure against the third virial term, 
which gives h ^ uj^ £^ / F^/'^ ^ f^^^'^- This is not shown 
in Fig. El 

In a good solvent, the polyelectrolyte brushes do not 
show any major qualitative difference from a theta sol- 
vent. A typical phase diagram is shown in Fig. [3 In this 
case, the first-order phase boundary occurs at higher val- 
ues of e and the critical point shifts to even smaller values 
of Pel- This is not surprising because in a good solvent, 
there is an additional outward pressure which the fluc- 
tuation pressure must overcome before the brush is col- 
lapsed. In addition, it is interesting to observe in Fig. |S1 
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FIG. 7: Typical phase diagram for good solvent {v = 17.0275, 
Z = 3) for the brush of figure|S] The solid line is a first order 
phase transition from the osmotic to collapsed brush regimes. 
For comparison, the transition for a brush in theta solvent is 
shown by the dotted line. Notice that in good solvent, the 
transition line occurs at higher values of e and the critical 
point occurs for smaller values of /3ej. This is also typical of 
brushes with large values of cu. As a further reference, the 
dashed line shows the curve of increasing charge fraction in 
for trivalent salt. For the case of monovalent counterions, the 
first-order transition occurs for smaller values of e than with 
trivalent counterions. 
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FIG. 8: Brush height vs. charge fraction in good solvent for 
monovalent counterions at high density {pa^ = 0.0625). 



where we have plotted the behavior of the brush height 
as a function of the charge fraction, that the brush height 
increases at low charge fraction, then decreases with in- 
creasing charge fraction, and eventually collapse. This 
nonmonotonic behavior is qualitatively similar to that 
observed by Bendejacq et al. [23| for monovalent salt in 
highly charged and dense cylindrical brushes. 

For an osmotic brush at the level of a scaling theory, 
the brush height is independent of the grafting density. In 
theta solvent, significant deviations from osmotic brush 
scaling can be seen even for monovalent counterions as 
well as a discontinuous collapse (see dotted curve in Fig. 
El). In good solvent, however, the brush height scales 
nearly linear in the grafting density, as shown in Fig. 
|51 This scaling can be understood from balancing of the 
second virial pressure against the fluctuation pressure: 
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FIG. 9: Brush height vs. density {N = 200, / = 0.3). 
In theta solvent (dotted curve), the brush height exhibits a 
first order collapse. In good solvent for monovalent (solid) 
and divalent (dashed) counterions, the brush height decreases 
initially, then increases. At large densities, the increase is 
linear. 



IV. CONCLUDING REMARKS 

In this paper, we have shown that expanding the 
counterion free energy systematically about Poisson- 
Boltzmann mean field theory for a single brush results 
in a collapse transition of the brush height with multiva- 
lent counterions. Though we find a first order transition 
similar to that of Csajka et al. [2^, we also find a critical 
point in the phase diagram and find that, generically, the 
decrease in brush height is continuous. 

The major difference in formulation between our the- 
ory and this scaling theory is that the fluctuations of 
counterions outside of the brush have been taken into 
account at the quadratic level. Fluctuations of counteri- 
ons from the interior to exterior region of the brush are 
likely to be responsible for the existence of the critical 
point in the phase diagram. 

Throughout this paper, we have assumed that all 
chains in the brush are equally stretched, and that the 
monomer density is constant throughout the brush. We 
have also employed the constant density approximation 
for the counterion density within the brush when expand- 
ing in the fluctuations. These approximations have been 
necessary in order to calculate the fluctuation pressure, 
though in reality the monomer density is known to be 
approximately parabolic 32]. It is possible that the dis- 
continuity of monomer density in our approximation is 
responsible for the flrst order transition region, and that 
a more exact theory may have a critical point at much 
smaller values of f3ei or no first order transition at all |33j . 

Neglected in our treatment is fluctuations of the poly- 
electrolyte chains themselves. These corrections to the 
free energy are likely to be unimportant if the fluctua- 
tions of uncondensed counterions do not couple strongly 
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to them. They may play a role, however, when condensed 
counterions are properly taken into account. 

We have also neglected the inhomogeneities of the 
monomer density in the brush, which induces inhomo- 
geneities in the mean field counterion density. These 
inhomogeneities are due, in part, to counterion conden- 
sation on the polyelectrolyte chains. These condensed 
counterions are not free to fluctuate and do not con- 
tribute to the osmotic pressure due to the counterion 
entropy. Further, the condensed counterions can cause 
additional attractions between the chains of the brush 
due to structual correlations beyond those calculated in 
this paper. 

In future work, we will study the problem of attractive 
interactions between two opposing brushes, and work to 
understand more precisely the role of condensed counte- 
rions in the brush. 
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APPENDIX A: CALCULATION OF 
FLUCTUATION ENERGY FOR A SINGLE 
BRUSH 

Within the constant counterion density approximation 
within the brush, the Fourier transform of G(x, x') in the 
xy-plane, G{z,z';q), satisfies 



(9^ 

-—^ + q^+A^e{z-h) + 



m^-h) 1 , 

^Ib5{z~z'), (A1) 



We acknowledge discussions with D. Bendejacq, V. 
Ponsinet, and in particular T.C. Lubensky and P. Pin- 



where A^ = 2 (1 — ^/A) / (^/i), and can be solved by stan- 
dard technique 34] . The result is 



G<(z,z;g) 



2a 



G>{z,z;q) = ^ 
2q 



1 



l + Mi{q)M2{q) 
1 1 ^i3/:(g)e-29(^-'«) 



2q 



1 



1 



q{z-h + \) 



for z < h. 



for z > h. 



^, , [1 - gA(l - qX)] {1 + + M2{q)e- 



2ah 



[1 -/-(?)]} 



(A2) 
(A3) 
(A4) 



[1 + qX{l + qX)] {1 + f+{q) + X2(g)e-2"^ [1 - /+(<z)]} ' 

where a'^ = q^ + A^, Mi{q) = e-^°'^ \ \+ull] ]' -^^l?) = (a - q)/{a + q), and /± = aX{l±qX)/ [l±qX{l±qX)] . 
Now, returning to the expression of the derivative of AjF with respect to h, Eq. H25|) . by noting that the derivative of 
dQ{z — h)/dh = —6{z — h) and dQ{z — h)/dh = +S(z — h), we find that Eq. H25|) consists of three parts: 

1 dAT 1 dATi 1 dAT2 1 dAT^ 



A dh A dh A dh A dh 



(A5) 



where A is the area of the plane and 

1 d(3ATi 
A dh 
1 d[3AT2 
A dh 
1 d!3AT3 
A dh 

After lengthy algebra, we have 

1 dAT A2 
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(2' 



= A'- 



A2 



°° dz d 

2i;^>'^'^'^'^dh 
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A dh h ' ^ ' -V /J ' x^hdh 
where the dimensionaless integrals Ti{h) are defined as 
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